Topological origins of a bi-parameter periodicity hub for the Rossler attractor 
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We explore the dynamical and topological characteristics of the Rossler system that lead to the 
existence of a periodicity hub and nested spiral in codimension-2 parameter space. We find that the 
hub shape is a consequence of conjugacy between the Rossler system and unimodal maps in the two- 
branch region of the parameter space. The nested spiral structure is a consequence of a topological 
feature of the Rossler system that has not been noted previously. We outline this mechanism and 
detail the spiral transition for the symbolic sequence of orbits up to period seven. 

PACS numbers: XXXXXXX 



Recent work has focused attention on the appear- 
ance of nested spirals and so called "shrimp" within the 
codimension-2 parameter space of a number of dissipa- 
tive systems. First noted by Bonatto and Gallas in 2008 
PQ , these hubs have been the focus of number of investi- 
gations since [2HE]. 

The occurrence of these hubs is easily seen in a map- 
ping of the global Lyapunov exponents in codimension-2 
parameter space (Figure [TJ. We focus on a hub found 
in the Rossler system. The hubs for this system have 
been studied for a number of different parameter values 
[2]. We chose to work with the parameter space of the 
Rossler system where b = 0.2 and (a, c) vary. It can be 
shown that the shape of this hub is robust with a varia- 
tion of 6, drifting rightwards and slightly downwards in 
the a, b parameter plane as b increases from 0.2 through 
2.0. 

The equations for the Rossler system [10] are, 

x = -y - z, y = x + ay, z = b + z(x-c) (1) 

While this system was first designed as a model for chaos, 
with no practical use in mind, numerous dynamical sys- 
tems in fields such as chemistry, electronics, biology, and 
lasers have been found to exhibit behavior conjugate to 
that of the Rossler system [TT] . 

In Figure [I] we draw an approximate best-fit curve 
linking the centers of so called "shrimp" structures along 
a line which corresponds to a transition from "spiral" (the 
return-map of the Rossler attractor has two branches) to 
"screw-like" structure in the Rossler attractor, which cor- 
responds to the branched-manifold describing the Rossler 
attractor obtaining a third branch and the return-map 
appearing bimodal with two critical points. We call this 
line the Topological Transition Line and write this partic- 
ular line as TTL\ to indicate the transition from needing 
a branched-manifold with two branches to one with three 
branches as we move in parameter space from left to right 
[12]. In agreement with previous estimates [6 , we find 
the center of the primary spiral hub to be numerically lo- 
cated at approximately (a, c) = (0.1798,10.3084) which 
corresponds to a homoclinic bifurcation point, the only 
such point along TTL\ pQ . 

We first wish to understand the smooth structure left 
of the TTL|, and we do so by assuming topological con- 
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FIG. 1: Lyapunov diagram for the Rossler system, close up 
of primary spiral hub. Green center dot (color online) indi- 
cates the center of the spiral hub. Dark blue line indicates a 
best-fit curve for the partition of the primary hub by the ap- 
pearance of "shrimp". Blue coloration reflects the intensity 
of the first Lyapunov exponent and corresponds to chaotic 
regions, red coloring reflects that of the negative of the sec- 
ond exponent when Ai = and corresponds to stable regions. 
Lyapunov exponents were calculated on a 1000x1000 grid us- 
ing the LESNLS software package [9] . 



jugacy between the Rossler system and unimodal maps 
for judiciously selected slices of parameter space. That 
is, while the return map for the Rossler attractor appears 
conjugate with unimodal maps left of TTL3, we seek a 
more specific conjugacy-that of the spectrum of stable 
periodic orbits. For this purpose, we will compare the 
Rossler system with the most widely studied unimodal 
map. The Logistic Map is described by the equation 
L(x n ,r) = x n+ i = rx n (1 — x n ). A method of analyti- 
cally detecting periodic windows in the general class of 
Logistic type maps was introduced by the seminal work of 
Metropolis et al. in 1973 [T3] , and rediscovered by Jensen 
and Myers in 1985 [El [15]. In order to locate windows of 
stability, we can study the iterations of the critical point. 
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Whenever one of these iterations of the critical point in- 
tersects the critical point, we are assured of the existence 
of a super-stable periodic orbit which corresponds to one 
of the windows of stability seen in the Logistic Map's bi- 
furcation diagram. These iterations of the critical point 
"scar" the bifurcation diagram as continuous "caustics" 
and stand out due to a higher concentration of "hits" 
around these points. A further advantage of the caustic 
method is that this method also describes the symbolic 
encoding of the super-stable orbits [13 . Similar caustics 
are seen in bifurcation diagrams taken for the Rossler 
system, but not necessarily with the same patterns as in 
the Logistic Map. 

The apparent conjugacy between the Rossler bifurca- 
tion diagram, and that of the Logistic Map, has not been 
commonly noted (see Letellier 1995 [16 for an exception). 
This is likely due to the fact that this similarity can only 
be found if one takes a slice through parameter space that 
is properly oriented with regards to periodicity conjugacy 
requirements between the maps. In order to judiciously 
select a slice in parameter space through which we can 
assume topological conjugacy with the spectrum of stable 
periodic orbits in unimodal dynamics, we note that the 
main conjugacy requirements of the Logistic Map will be 

• Unimodal structure. 

• An initial period- 1 orbit that undergoes a period dou- 
bling cascade and then transitions to chaos upon reach- 
ing the Feigenbaum point. 

• The attractor evolves until it makes contact with a fixed 
point. 

The Rossler system matches these requirements in a spe- 
cific region, 

• Unimodal structure: due to the highly dissipative na- 
ture of the Rossler system, its return map has a uni- 
modal nature in the regions left of TTL\. 

• The Rossler system undergoes an Andronov-Hopf Bi- 
furcation, which leads to a period- 1 limit cycle that 
undergoes a period-doubling cascade as we get closer 
to TTL\ in parameter space. 

• The Rossler system has a homoclinic bifurcation curve 
(a thin parabola shape), the tip of which lies on TTL\ 
PP- The attractor makes contact with its fixed point in 
the form of a homoclinic orbit at this point. 

We treat the Rossler system as if it were fully topolog- 
ically conjugate to the Logistic Map for any line in the 
parameter space which begins at the Andronov-Hopf Bi- 
furcation and terminates at the homoclinic point on the 
TTL\ where the Rossler attractor first makes contact 
with its central fixed point. This results in the shape 
of the hub left of TTL\ tracing out the shape of the 
Andronov-Hopf bifurcation curve. 

In Figure [2j we draw an example of two such slices of 
parameter space that will fulfill these conjugacy require- 
ments, LI and L2. We find that their bifurcation diagram 
is qualitatively similar to that of the Logistic map, and 
furthermore, we tested this assumption by numerically 
recreating the first seven caustics in these bifurcation di- 
agrams. At least up to period seven and likely so for 



higher periods, the order of periodic orbits in the Rossler 
system as we move along these lines is conjugate with 
that of the Logistic Map. To further test these assump- 
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FIG. 2: Parameter space diagram with Andronov-Hopf curve 
calculated using PyCONT/AUTO [13 [18]. Lines LI and 
L2 originate from the period- 1 limit cycle created after the 
Andronov-Hopf bifurcation of the Rossler attractor and termi- 
nate at a Homoclinic bifurcation point on the curve marking 
the topological transition from two to three branches in the 
branched manifold. The order of stable periodic orbits seen 
along these lines are conjugate to those seen in the Logistic 
map up to at least period seven and likely for much higher 
periods. We exclude the continuation of the homoclinic bi- 
furcation curve into the three-branch region for simplicity. 




FIG. 3: Periodicity calculations for the codimensional-2 pa- 
rameter space of the Rossler attractor, for b = 0.2, a G 
(0.1,0.22), c G (5.0,25.0). Darkest blue corresponds to di- 
vergent or chaotic behavior. 

tions, we partitioned parameter space around the spiral 
hub into a 5000 by 5000 grid and applied a periodicity 
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detection algorithm for each point. We color code the 
periods, up to period 11, producing an extensive period- 
icity graph in Figure [3] which shows that the periodicity 
is structured as predicted. 

Finally, we can see that by continuity, this ordering 
of periods also determines the order of periods of the 
shrimp structures along the TTL\. At the center of these 
shrimps, one finds a degenerate super-stable state which 
corresponds to the fact that at these points two intersect- 
ing parabola-shaped regions of super-stable periodicity 
contain both critical points simultaneously [15] . Care- 
ful consideration of the conditions at the intersection of 
the topological regions will help us understand why the 
shrimp line up along the transition curve. These stable 
regions represent iterations of at least one critical point 
of the return-map. The first two iterations of the critical 
point mark the borders of the range of possible itera- 
tions. As we pass through TTL\ in phase space from 
left to right, the Rossler's state-space return map grows 
a third branch. At TTL|, the edge of the second branch 
is precisely the second critical point of the bimodal map 
in the branch-3 region. 




FIG. 4: Anatomy of a generic "mutant shrimp" along the 
largest spiral. The symbolic string describing these orbits 
will be dominated by zeros. The center of the shrimp struc- 
ture corresponds to a doubly-superstable orbit which includes 
both critical points. Each parabola corresponds to a varia- 
tion of this orbit in which only one of the critical points is 
on the orbit. Due to the topological mechanism responsible 
for the periodicity transition, this shrimp has one tail which 
originates as a lower period and transitions to higher period 
before crossing TTL\. 

The shrimp along the TTL\ and above the focal point 
behave as we would expect in the bimodal regime [15] . 
One of the "tails" of each of the shrimp points downward 
back towards the hub. This "tail" joins a shrimp below 
the focal point which has a different period. We call 
these "mutant shrimp" and have detected the periodicity 



transition to happen very close to TTL\ (Figures p] and 

%. 

Directly related to these structures is the question of 
the source of the nested periodicity spirals. These spirals 
are due to a previously unnoticed topological feature of 
the Rossler system. The branched-manifold describing 




FIG. 5: As we trace the periodicity spiral around in param- 
eter space, clockwise from lower to higher period, the third 
branch of branched-manifold in state space is reinjected on 
the plane of the attractor at an angle that varies clockwise 
around the fixed-point, merging with the second branch as it 
crosses the TTL\ from right to left. 

the Rossler system has been typically assumed to have 
a standard shape in which all branches are redeposited 
onto the main plane of the attractor at approximately 
the same position. We find that this is not the case for a 
small range of values in codimension-2 parameter space. 
In particular, as we move around one of the spiral fea- 
tures in the hub in parameter space, we find that in state 
space, the position at which the third branch is deposited 
on the main plane of the attractor rotates around the 
fixed point. Figure [5] shows this clearly for the spiral 
transition from period three to period four. As we follow 
a stable periodic orbit clockwise around the central spi- 
ral hub from lower periodicity to higher periodicity, the 
angle of deposit also rotates clockwise around the central 
dynamical fixed point in state space. Like a long chain 
being deposited in a circular way on a flat surface, the 
number of turns of the chain on the flat surface increases 
as we follow the periodicity regions as they spiral in to- 
wards the center of the spiral hub. As we cross TTL\ 
following these spirals, the third branch "collapses" back 
onto the second branch, and the cycle starts over. This 
results in the transition between lower to higher period 
happening very rapidly near the shrimp along the bottom 
section of TTL\ below the homoclinic point. 

Following these spirals inwards, we find that the new 
points on an orbit are deposited on the inner-most region 
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of the attractor, winding the attractor closer and closer 
to its fixed point until it makes full contact at the ho- 
moclinic bifurcation. Knowing this new mechanism, we 
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FIG. 6: Symbol sequences in the nested spiral hub up to 
period seven. The orbits listed in each column are in the order 
that they appear in unimodal mappings [12]. Green arrows 
indicate orderings which have been matched both by looking 
at the structure in state space, and by matching the symbol 
sequence of the corresponding transitions in parameter space. 
The transition CD0111 CD00111 has only been confirmed 
by visually inspecting the attractor. Orbits in bold originate 
from isolated shrimp off of TTL\ within the branch-3 side. 
The attached gray boxes indicate the approximate position of 
the corresponding mutant shrimp in parameter space. 



Thus both connected periods p and p+1 will begin with 
symbol C. For example, we know that the period-four 
orbit C100 must transition into an orbit CTX00 after 
passing through the sequence C200 and returning to the 
TTL\ since the third branch lays the new point on the 
orbit, where X is either or 1. Since we know that the 
new branch is deposited on the inner region of the main 
plane of the attractor, this new symbol will be a 0, and so 
we have the transition from C100 C200 C1000. We 
have traced the trains for a number of orbits (Figure [6| 
which follow this general algorithm, though we note that 
higher periods will have a more complicated algorithm. 
We also find that the orbit CDO is connected with an 
orbit C2 inside the branch-3 region, and CD01 has been 
found to be connected to an orbit C21 also from this re- 
gion. Finally, we note that CD01 is the period-double of 
the orbit CI which itself is connected to a period-one or- 
bit in the branch-3 region which transitions to the orbit 
CI by the same mechanism outlined in this report. The 
shape of this period- 1 orbit is outlined in the upper left 
of Figure [6] as it nears transition. 

Orbits that can not originate within a spiral, because 
the would-be preceding orbit does not exist, instead orig- 
inate from one of the isolated shrimp in the branch-3 side 
of parameter space, though they themselves are then the 
beginning point of a new spiral, resulting in a nested set 
of spirals. 



can predict the symbolic sequence of orbits along a spi- 
ral. We must have that the critical point associated with 
the unimodal map on the branch-2 side of TTX§, which 
we label C, is always part of the symbolic sequence along 
this spiral. This is due to the fact that the spiral region is 
continuous and crosses over the branch-2 side where the 
bimodal critical point, which we label D, is excluded. 
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